SPECIFICATION 



TO ALL WHOM IT MAY CONCERN: 

Be it known that I, Michael Kallay, a citizen of the 
United States, residing at 9815 Vineyard Crest, Bellevue, 
Washington 98004, have invented a certain new and useful 
METHOD AND SYSTEM FOR COMPUTING THE MASS PROPERTIES OF 
POLYGONS AND POLYHEDRA of which the following is a 
specification . 



METHOD AND SYSTEM FOR COMPUTING THE MASS PROPERTIES 
OF POLYGONS AND POLYHEDRA 



5 FIELD OF THE INVENTION 

The invention relates generally to computer systems, and 
more particularly to computer graphics. 



BACKGROUND OF THE INVENTION 

10 The center of mass and the moment of inertia tensor are 

known as the mass properties of a rigid body. Engineers use 
mass properties them to analyze the strength of mechanical 
structures, and to compute the dynamics of moving bodies. 
Similarly, mass properties can be used in computer 

15 simulations. As physics-based simulation is becoming 

important in movies and computer games, the mass properties 
are gaining importance in both two-dimensional and three- 
dimensional computer graphics. For example, to realistically 
simulate the motion of a moving body following a collision, 

20 the mass properties may be used to determine where the object 
should appear in subsequent frames, e.g., the moment of 
inertia of a body is a primary factor in determining the speed 
of the rotation of the body following a collision. 

While computing the center of mass is straightforward, 

25 computing the moment of inertia is much more elaborate. 

Standard mathematical references contain formulae for some 



special shapes (e.g. rectangles, cylinders) in terms of their 
dimensions. For polyhedrons in general, an algorithm for 
computing the integrals on a polyhedron has been presented in 
the reference by Brian Mirtich, Fast and Accurate Computation 
5 of Polyhedral Mass Properties, Journal of Graphics Tools, 
volume 1, number 2, 1996 (the Mirtich reference) . This 
algorithm relies on a scheme proposed in Sheung-ling Lien and 
James T. Kajiya, A symbolic method for calculating the 
Integral properties of arbitrary nonconvex polyhedra, IEEE 

10 Computer Graphics and Applications, 4(10), pp 35-41, 1984, 

which describes converting a volume integral on a polyhedron 
to a sum of surface integrals on its facets, using the 
Divergence theorem. That scheme is used in the aforementioned 
Mirtich reference, to convert the resulting facets integrals 

15 to line integrals on their edges, using Green's theorem. 

As the images rendered by computer graphics become more 
realistic, significant numbers of computations at high frame 
rates are required, however large numbers of overly complex 
calculations (e.g., on hundreds or even thousands of polygons) 

20 may not be able to be performed fast enough on commercially 
available computing devices such as gaming devices. Further, 
reducing the complexity so that developers can better use 
quadratic homogeneous polynomials (such as the entries of the 
moment of inertia tensor) is also desirable. What is needed 



is a more straightforward way to compute the mass properties 
of polygons and polyhedra. 

SUMMBMRY OF THE IMVENTION 

5 Briefly, the present invention provides a method and 

system in which formulae express integrals on polygons and 
polyhedra in terms of vertex coordinates. As a result, the 
moment of inertia of a polygon can be computed without 
integrating, e»g., for a triangle by using multiplication of 
10 one-twelfth the area times the sum of the dot products of (or 
based on) the vertices. 

In one implementation, the entries of the moment of 
inertia tensor and the moment of inertia of a body about a 
given axis through the origin are volume integrals of 
15 homogeneous quadratic polynomials. 

If r is a triangle with vertices ^,,^2*^*3 then the 
following two formulae apply: 

(3) ff/(X-C) = -^t(/(/^)-g(P,.„/^)), 

/ lo /^i 

20 and 

(4) JJ/(Jr-c)=^t/(i^.-c), 
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where f{X) is any homogeneous quadratic function of the 
coordinates, giXy.X^) is its polar form, a{T) is the area of T, 
and C is the centroid of T. Given these formulae, it is seen 
that the moment of inertia (or like quadratic homogeneous 
5 polynomials) can be determined without integration. For 

example, with a triangle of vertices Px^Pi.P^ with an origin at 

its centroid, the moment of inertia is ^^(P, -/j + A '^3 +^3 "^i) • 

1. 2 

Thus, simple multiplication and division (applied to the sums 
of the dot products of the vertices) provides the solution. 

10 Moreover, as is known, particularly in computer graphics 

technology, any polygon may be divided into a sum of 
triangles. Thus, the above present invention may be applied 
to any polygon, including one with holes. A simple 
triangulation process may establish the triangles as is known 

15 in the art. However, if B is non-convex or does not contain 
the origin, then some of these triangles may overlap, while 
some may spill out of B . Nevertheless, equality is preserved 
through the use of a signed sum of the triangles, where the 
sign is determined by the orientation (clockwise or 

20 counterclockwise) of the path Any overlapping and/or 

spilling out is handled in a straightforward way by replacing 
the area of each triangle with one-half of the determinant 
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relative to an origin, Y^et(P^.,,P;) , which inherently handles the 

proper sign. Taking the sum of the signed triangles that make 
up any polygon provides the solution • As a result, 
integration can be avoided regardless of the shape of the 
5 polygon. 

Various aspects of the present invention may be applied 
to three dimensional polyhedra via integrating over a cone. A 
general polyhedron is a '^signed union'^ of cones extended from 
an origin to its facets, where the sign is determined by the 

10 side of the plane of the facet that contains the origin, 

positive if the origin lies on the interior side. When the 
body is a simplicial polyhedron, i.e., when all its facets are 
triangles, which is common in computer graphics, if every 
facet of a body 5 is a triangle with vertices/^.,, , and 

15 /is a homogeneous quadratic polynomial then: 



where g is the polar form of / . 

Thus, a general polyhedron may be made up of the signed 
20 sum of the cones / pyramids. Thus, for a general polyhedron, 
in which the body 5 is a polyhedron with facets F,,,,F„, and / 



(10) 
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10 



is a homogeneous quadratic polynomial whose polar form is g 
then JJ//W^ = -}2W ^^/M^)(/"(C,) + 

TVZdet(^,:^.iv,)(/-(^;;:;)+/(^)+g(^,^)) 

where C, is the centroid and aiF^) is the area of , and 
^denotes /^,;-C, - 

Note that the centroid and area of F/ can be obtained 



from: 



(12) fl(F,) = i|;'det(P,,_, -P,„P,j -P,„N,) 



where is a unit vector normal the facet F., . 

Any mass properties formula can now be generated by 
defining the appropriate function / and multiplying by 
density. For example, \ii.tb. f {X) = X ■ X , and its polar form 
15 g{X^,X2) = X^■X2 we have from formula (4) that the moment of 
inertia of a triangle with vertices P^,P^,P^, mass m and 
centroid C, about an axis perpendicular to the plane of T at C 



is: 
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^t(^,-C).(i^.-C), 

or: 



5 In keeping with the present invention, it follows from 

formula (5) that the moment of inertia of a triangle with 
vertices P^^P2,P^ and mass m about an axis iO,V) is: 

1 {{P^ X V) - (i^. X F) + (P,_, X V) . (P, X n) . 
1=1 



The integrands for the entries of moment of inertia 



tensor 


are: 






For I„ : 


/W = 


j'^+z^ and the polar form is^(X,,^2) 




For/_^ : 


/W = 


z^+x^ and the polar form is g(X^,X2) 


= ZjZ2 + ^1^2 • 


For/^ : 




x^+y^ and the polar form isg(X,,X2) 




For/^: 




-xy and the polar form isg(X^yX2)-- 




For/^: 


/w= 


-yz and the polar form ±s g{X^^X2)-- 




For/„ : 




-zr and the polar form isg(XpX2) = - 


^(2,X2+X,Z2) . 



The present invention can selectively combine the 
formulae to produce a desired result for any polyhedron. For 
20 example, after initializing the resulting integrals (which 
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will be summed) to zero, for each facet with outward normal 
Ni, if the facet is a triangle (i^,,,/^,2.^3) / the area of the 
triangle is computed from the vertices, using the above 
described formula a(T) = i|det(i»2 -/^ ,,P,^ , while centroid is 
5 computed as C, = +i^,2 +^.3) • 

With the centroid and area known, and applied to solve 



the 



formula jjfiX -Q = ^^"^f(P. -C) , the appropriate facet 



integrals about centroid C, are easily computed. 

If the facet is not a triangle, the centroid is computed 



10 



as C,=^^Xdet(^.>-.-^..'^j-^..'^/ Then, the 

appropriate facet integrals about C, are computed using 
jjnX)da = jkl;det(/^.-„i^.)(/(P,.,) + /(/^.) + ^(i^..,.P,)) . 

With the integrals about the centroid of the facet C,., 
whether computed for a triangle or otherwise, the facet 
15 integrals about the origin are computed using 

l\f(X)da= jjfiX-C)da + a(B)f(C). The signed volume integrals over 

B B 

the facet's cone are computed using jljf{X)(iv^\h(F)jjf{X)da , 

B F 

with iV. for the signed height. The computed integrals are 
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added to the accumulative resulting integrals, until no facets 
remain and the accumulative resulting integrals are known. 

In one example implementation, the formulae may be 
incorporated in a graphics computation process. Because 
5 contemporary computer graphics use large quantities of 

triangles, thus the present invention can provide substantial 
benefits when the mass properties of an object or multiple 
objects need to be recomputed. 

Other advantages will become apparent from the following 
10 detailed description when taken in conjunction with the 
drawings, in which: 



BRIEF DESCRIPTION OF THE DRaWINGS 

FIGURE 1 is a block diagram generally representing a 
15 computer system into which the present invention may be 
incorporated; 

FIG. 2 is a representation of a triangle divided into 
four triangles by connecting the midpoints of its edges in 
accordance with an aspect of the present invention; 
20 FIG. 3 is a representation of integrating over a 

polyhedron, in accordance with an aspect of the present 
invention; 
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FIGS- 4 and 5 comprise a flow diagram generally 
representing the operations for integrating over a polyhedron, 
in accordance with an aspect of the present invention; and 

FIG. 6 is a block diagram generally representing an 
5 example graphics processing architecture into which the 
present invention may be incorporated. 



DETAILED DESCRIPTION 

EXEMPLARY OPERATING ENVIRONMENT 

10 FIGURE 1 illustrates an example of a suitable computing 

system environment 100 on which the invention may be 
implemented. The computing system environment 100 is only one 
example of a suitable computing environment and is not 
intended to suggest any limitation as to the scope of use or 

15 functionality of the invention. Neither should the computing 
environment 100 be interpreted as having any dependency or 
requirement relating to any one or combination of components 
illustrated in the exemplary operating environment 100. 

The invention is operational with numerous other general 

20 purpose or special purpose computing system environments or 
configurations. Examples of well known computing systems, 
environments, and/or configurations that may be suitable for 
use with the invention include, but are not limited to: 
personal computers, server computers, hand-held or laptop 
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devices, tablet devices, multiprocessor systems, 
microprocessor-based systems, set top boxes, programmable 
consumer electronics, network PCs, minicomputers, mainframe 
computers, distributed computing environments that include any 
5 of the above systems or devices, and the like. 

The invention may be described in the general context of 
computer-executable instructions, such as program modules, 
being executed by a computer. Generally, program modules 
include routines, programs, objects, components, data 

10 structures, and so forth, which perform particular tasks or 
implement particular abstract data types. The invention may 
also be practiced in distributed computing environments where 
tasks are performed by remote processing devices that are 
linked through a communications network. In a distributed 

15 computing environment, program modules may be located in local 
and/or remote computer storage media including memory storage 
devices . 

With reference to FIG. 1, an exemplary system for 
implementing the invention includes a general purpose 
20 computing device in the form of a computer 110. Components of 
the computer 110 may include, but are not limited to, a 
processing unit 120, a system memory 130, and a system bus 121 
that couples various system components including the system 
memory to the processing unit 120. The system bus 121 may be 
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any of several types of bus structures including a memory bus 
or memory controller, a peripheral bus, and a local bus using 
any of a variety of bus architectures. By way of example, and 
not limitation, such architectures include Industry Standard 
5 Architecture (ISA) bus. Micro Channel Architecture (MCA) bus. 
Enhanced ISA (EISA) bus. Video Electronics Standards 
Association (VESA) local bus, and Peripheral Component 
Interconnect (PCI) bus also known as Mezzanine bus. 

The computer 110 typically includes a variety of 

10 computer-readable media. Computer-readable media can be any 
available media that can be accessed by the computer 110 and 
includes both volatile and nonvolatile media, and removable 
and non-removable media. By way of example, and not 
limitation, computer-readable media may comprise computer 

15 storage media and communication media. Computer storage media 
includes volatile and nonvolatile, removable and non-removable 
media implemented in any method or technology for storage of 
information such as computer-readable instructions, data 
structures, program modules or other data. Computer storage 

20 media includes, but is not limited to, RAM, ROM, EEPROM, flash 
memory or other memory technology, CD-ROM, digital versatile 
disks (DVD) or other optical disk storage, magnetic cassettes, 
magnetic tape, magnetic disk storage or other magnetic storage 
devices, or any. other medium which can be used to store the 
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desired information and which can accessed by the computer 
110. Communication media typically embodies computer-readable 
instructions, data structures, program modules or other data 
in a modulated data signal such as a carrier wave or other 
5 transport mechanism and includes any information delivery 
media. The term ^^modulated data signal'' means a signal that 
has one or more of its characteristics set or changed in such 
a manner as to encode information in the signal. By way of 
example, and not limitation, communication media includes 
10 wired media such as a wired network or direct-wired 

connection, and wireless media such as acoustic, RF, infrared 
and other wireless media. Combinations of the any of the 
above should also be included within the scope of computer- 
readable media. 

15 The system memory 130 includes computer storage media in 

the form of volatile and/or nonvolatile memory such as read 
only memory (ROM) 131 and random access memory (RAM) 132. A 
basic input/output system 133 (BIOS), containing the basic 
routines that help to transfer information between elements 

20 within computer 110, such as during start-up, is typically 
stored in ROM 131. RAM 132 typically contains data and/or 
program modules that are immediately accessible to and/or 
presently being operated on by processing unit 120. By way of 
example, and not limitation, FIG. 1 illustrates operating 
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system 134^ application programs 135, other program modules 
136 and program data 137. 

The computer 110 may also include other removable/non- 
removable, volatile/nonvolatile computer storage media. By 
5 way of example only, FIG. 1 illustrates a hard disk drive 141 
that reads from or writes to non-removable, nonvolatile 
magnetic media, a magnetic disk drive 151 that reads from or 
writes to a removable, nonvolatile magnetic disk 152, and an 
optical disk drive 155 that reads from or writes to a 

10 removable, nonvolatile optical disk 156 such as a CD ROM or 
other optical media. Other removable/non-removable, 
volatile/nonvolatile computer storage media that can be used 
in the exemplary operating environment include, but are not 
limited to, magnetic tape cassettes, flash memory cards, 

15 digital versatile disks, digital video tape, solid state RAM, 
solid state ROM, and the like. The hard disk drive 141 is 
typically connected to the system bus 121 through a non- 
removable memory interface such as interface 140, and magnetic 
disk drive 151 and optical disk drive 155 are typically 

20 connected to the system bus 121 by a removable memory 
interface, such as interface 150. 

The drives and their associated computer storage media, 
discussed above and illustrated in FIG. 1, provide storage of 
computer-readable instructions, data structures, program 
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modules and other data for the computer 110. In FIG. 1, for 
example, hard disk drive 141 is illustrated as storing 
operating system 144, application programs 145, other program 
modules 146 and program data 147. Note that these components 
5 can either be the same as or different from operating system 
134, application programs 135, other program modules 136, and 
program data 137. Operating system 144, application programs 
145, other program modules 14 6, and program data 147 are given 
different numbers herein to illustrate that, at a minimum, 

10 they are different copies. A user may enter commands and 

information into the computer 20 through input devices such as 
a tablet, or electronic digitizer, 164, a microphone 163, a 
keyboard 162 and pointing device 161, commonly referred to as 
mouse, trackball or touch pad. Other input devices not shown 

15 in FIG. 1 may include a joystick, game pad, satellite dish, 

scanner, or the like. These and other input devices are often 
connected to the processing unit 120 through a user input 
interface 160 that is coupled to the system bus, but may be 
connected by other interface and bus structures, such as a 

20 parallel port, game port or a universal serial bus (USB) . A 
monitor 191 or other type of display device is also connected 
to the system bus 121 via an interface, such as a video 
interface 190. The monitor 191 may also be integrated with a 
touch-screen panel or the like. Note that the monitor and/or 
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touch screen panel can be physically coupled to a housing in 
which the computing device 110 is incorporated, such as in a 
tablet-type personal computer • In addition, computers such as 
the computing device 110 may also include other peripheral 
5 output devices such as speakers 195 and printer 196, which may 
be connected through an output peripheral interface 194 or the 
like. 

The computer 110 may operate in a networked environment 
using logical connections to one or more remote computers, 

10 such as a remote computer 180. The remote computer 180 may be 
a personal computer, a server, a router, a network PC, a peer 
device or other common network node, and typically includes 
many or all of the elements described above relative to the 
computer 110, although only a memory storage device 181 has 

15 been illustrated in FIG. 1. The logical connections depicted 
in FIG. 1 include a local area network (LAN) 171 and a wide 
area network (WAN) 173, but may also include other networks. 
Such networking environments are commonplace in offices, 
enterprise-wide computer networks, intranets and the Internet. 

20 When used in a LAN networking environment, the computer 110 is 
connected to the LAN 171 through a network interface or 
adapter 170. When used in a WAN networking environment, the 
computer 110 typically includes a modem 172 or other means for 
establishing communications over the WAN 173, such as the 
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Internet. The modem 172, which may be internal or external, 
may be connected to the system bus 121 via the user input 
interface 160 or other appropriate mechanism. In a networked 
environment, program modules depicted relative to the computer 
5 110, or portions thereof, may be stored in the remote memory 
storage device. By way of example, and not limitation, FIG. 1 
illustrates remote application programs 185 as residing on 
memory device 181. It will be appreciated that the network 
connections shown are exemplary and other means of 
10 establishing a communications link between the computers may 
be used. 

COMPUTING THE MASS PROPERTIES OF POLYGONS AND 
POLYHEDRA 

The present invention is generally directed towards a 
method and system for expressing integrals of quadratic 
homogeneous polynomials (where ( f{tX) ^ f{X)) on polygons and 
polyhedra in terms of vertex coordinates. For example, the 
entries of the moment of inertia tensor are integrals of 
quadratic homogeneous polynomials, and thus the present 
invention is particularly beneficial in determining and 
expressing the entries of the moment of inertia tensor in 
terms of simple vertex coordinates. As can be readily 
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15 



20 



25 



appreciated, expressing complex mathematical formulae in the 
form of simple vertex coordinates facilitates efficient 
processing in many environments, such as computer graphics 
simulation in which many computations may have to be done 
5 within a relative short time frame. Notwithstanding, the 
present invention is not limited to computing the moment of 
inertia, nor to computer graphics processing, but rather 
provides many benefits in contemporary computing. 

The center of mass and the moment of inertia tensor are 
10 known as the mass properties of a rigid body. The center of 

mass of B is C jj^Xpdv , where/? is density and v is volume; 

B 

this is a vector notation, andX stands for {x^y.z) . The moment 
of inertia of 5 about an axis is the integral where r 

B 

is distance from the axis. Writing dm foxpdv, the entries of 
15 the moment of inertia tensor are: 

B B B 

ly^-l^' Iyy = W^^'dm , I^=-\\\yzdm, 

B B 
B 

20 
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When the body B is of uniform density, the density {p) 
is a constant that may be factored out, and the remaining 
integrals are purely geometric properties of B . The center of 
s is then called the centroid of B, which is C^\^Xdv. For 



mas 



5 the other properties, the density is assumed to be 1 (p=l) for 
computing the integrals. When 5 is a polyhedron or a planar 
polygon, closed formulae for these integrals are derived, in 
accordance with an aspect of the present invention. 

A planar polygonal plate 5 is a body that is cut out of a 
10 planar plate of thickness t in the shape of some region R . If 
t is sufficiently small, then the integrand changes very 
little (or not at all in some cases) along a direction 
perpendicular to the plate, and then the volume 

integral JJJ/(A")i/v is approximated by the area integral 



15 t\\f{X)da 



R 

In accordance with the present invention, a more 
geometric approach is provided in which the formulae for 
polygons are derived with elementary geometry, and slab 
integration is used to obtain the integral on a polyhedron. 
20 Note that a general framework for computing the integral of 
any homogeneous quadratic polynomial over polygons and 
polyhedra is presented herein, employing the concept of polar 
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forms, wherein a polar form for such a polynomial /(JT) is a 
bilinear symmetric function of two (vector) variables g(X^,X2) 
for which g(X,X) = f(X). The resulting formulae for the mass 
properties of triangles and polygons are more concise than 
those of those of the prior art, and thus are easier to 
implement, less error prone, and likely to lead to faster 
and/or more accurate algorithms. 

As used herein, O stands for the origin. For a body B, 
a real number t and a pointP, /^denotes the point-set 
{tXiXsB}, and B + P denotes the set {X-^PiXsB}. The 
''polygon(Pi,„P„)" stands for a polygon having vertices Pj„,P„ , 
with the convention Po=J^«- Similarly, ''the axis 
represents an axis through the point P that points in the 
direction of the unit vector V . A facet is a simple polygon 
with a consistently oriented sequence of vertices i^. ^^.^ , 
with the convention /^•o=^i,«(o- Each facet F. is assumed to 
possess a well-defined, out -pointing normal unit vector JV,.. 
The vertices are assumed to be oriented clockwise when looking 
in the direction ofiST. (i.e., from the interior of B), 

In one implementation, the entries of the moment of 
inertia tensor and the moment of inertia of a body about a 
given axis through the origin are volume integrals of 
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homogeneous quadratic polynomials. As described above, the 
polar form for such a polynomial f{X) is a bilinear symmetric 
function of two (vector) variables g{X^,X^) for which 
g{X,X)- f{X) . Although not necessary to the present invention, 
the existence and uniqueness of a polar form has been proved, 
and g can be referred to as the polar form of / . The 
linearity and symmetry of g imply the following formula: 

(1) f{P±Q) = fiP)^f{Q)±2g{P,Q) . 

The following is a well known generalization of the so 
called Parallel Axis Theorem: 

Lemma: If Cis the centroid otB , then: 
( 2 ) \\\KX)dv = \\\f{X - C)dv + v{B)f(Cr 

B B 



The planar version of formula (2) is another well known 
formula : 

(2a) \\f{X)da^\\f{X-C)da + a{B)f{C), 

B B 

where a{B) is the area of B . 
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For a triangle, the centroid of a triangle T with 



vertices/^, is simply {^Pf . As described below, the triangle 



also has relatively straightforward formulae for its mass 
properties. 

More particularly, consider the following theorem in 
which / is a homogeneous quadratic polynomial, and g is its 
polar form. If T is a triangle with vertices PpPj^^s then the 
following two formulae apply: 



where a{T) is the area of T, and C is the centroid of T. 

Proof: As represented in FIG. 2, T may be divided into 
four triangles T^^T^^T2,T^ by connecting the midpoints of its 
edges, giving four similar triangles, but having half-length 
edges. An integral over the large triangle would be the same 
as four times an integral over the small triangle, plus a term 




and 



(4) 
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3 

for the origin shift for each. Since r = (j7]. and the interiors 
of 7] are disjoint, then it follows that: 
jjfiX -C)da = j;^ \\f{X - C)da , 

and by formula (2) above, that: 



\\f{X-Oda = X J//(X-C,)flfa + fl(7;.)/(C, -C) , 

T i«=0 t; 



where C,- is the centroid of 7)- . However, these sub- 
10 triangles are half-size copies of T with respect to their 

sides, hence fl(7)) = {a(r). Moreover, as an area integral of a 

quadratic homogeneous polynomial, JJ/(Jf-C^)d!af = ^ J|/(jr-C)rffl • 

Ti T 

This is true for / = 1,2,3, where 7]. is a straight copy of ~T . It 
is somewhat not as apparent for T^, which is a copy of , 
15 but it is true since/is a quadratic homogeneous polynomial, 
hence /(~^) = /(X). Thus 

\\f{X ^ C)da = 1 \\f{X - C)da + \ a(T)f^ /(C, - C) . 
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This is an equation in jjf(X-C)da whose solution is; 

T 

|J/(X-CMa = ia(r)2/(C, -C) . 



Because Q - C = 0 , /(Q - C) = 0 . For / = 1,2,3 : 
5 C, =i(P, +!(/>_, + +i>.,)) = i(4/> +/^.., . 

Further, 

-C = i(4P, +/^,,)-K/^ +/^.i) = i(2/^ -/^.,) 

hence 

10 /(C,-C) = ^/(2/^-i^._,-i^.,,). 

Applying formula (1) and summing up, leaves: 



3 



x;/(c, -c) = 3^2:(6/(/^.)-6g(/^..,./:.))= \±{f(p,)-g(p,.„p,)), 



which proves formula (3) . Moreover 

C, =i(3/^. +(/^. +P,^,)) = i(3i^. +3C) = i(/^ +C) , 

15 hence 

/(C, - C) = /(I (i^. + C) - C) = i fiP, - C) , 
which proves formula (4). 

Given these formulae, it is seen that the moment of 
20 inertia (or like quadratic homogeneous polynomials) can be 

determined without integration. For example, with a triangle 
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of vertices P,,P2 7 with an origin at its center, the moment of 



and division (applied to the sums of the dot products of the 
vertices) provides the solution. 

Moreover, as is known, particularly in computer graphics 
technology, any polygon may be divided into a sum of 
triangles. Thus, the above present invention may be applied 
to any polygon, including one with holes. 

As another theorem, let / be a homogeneous quadratic 
polynomial and let g be its polar form. If T is a triangle 
with vertices P^,P2,P^ then: 



where a(T) is the area of T. 



Proof: Combining formula (2) and formula (3) results in: 



inertia is 



area 



(Pi'Pi+Pi'P^-^Ps'Pi)' Thus, simple multiplication 




jlnX)da = jinx ^ C)da + a(T)fiC) = 



3 3 



^a(7')5]/(/^)-g(/^..„/>) + ia(D/(£i^.) = 



/-I 1=1 



iLa(r)(^/(P,)-g(/'M,^.) + 2t(/(/^) + 2s:(/>..,/',)) 




3 



ia(r)2((/^)+g(/^-„/^)). 
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Turning to a consideration of integrating over a polygon, 
another theorem states: Let £ be a planar region in the(jc,>^) 
plane bounded by a simple counterclockwise polygon (P, P^) . 
Then, for any homogeneous quadratic polynomial / and its polar 
form g, another formula is: 

( 6) jjfiX)da = iL2det(/^._„i^)(/-(/>_,) + /(/^.) + g(i^_„i^.)) . 

B i=l 

Proof: If 5 is convex and it contains the origin O, 
then it is the disjoint union of triangles 7) with vertices 

Since f(O) = g{X,O) = 0 , it follows from formula (5) that: 
JJ/W = J//W = ii:a(7;.)(/(P,_,) + /(P,) + ^(/>.,,i^)) . 

B '=1 7) '=1 

A simple triangulation process may establish the 
triangles as is known in the art. However, if B is non-convex 
does not contain the origin, then some of these triangles may 
overlap, while some may spill out of B. Nevertheless, 
equality is preserved through the use of a signed sum of the 
triangles, where the sign is determined by the orientation 
(clockwise or counterclockwise) of the path . 

This overlapping and spilling out is handled in a 
straightforward way by replacing a{ri) with one-half of the 



determinant relative to an origin, ^det(P^_,,/;.) , which inherently 
handles the proper sign. Taking the sum of the signed 
triangles that make up any polygon provides the solution. 
The same argument applies to a planar polygon in any 
5 plane if replacing O with some point in the plane. For 

example, suppose the body B lies in a plane perpendicular to 
the unit vector N , and |g is a point in that plane. Denoting 

X^X-Q, the signed area of the triangle is 

ydet(/5_,,i^._j,A'^), hence for a clockwise orientation when viewed in 
10 the direction of N : 

( 7 ) \\f{X - Q)da = i^JdetC^,:^, A^)(/(^) + /(^) + giPT^JP,)) . 

B M 



In sum, it is shown that integration can be avoided 
regardless of the shape of the polygon. More complex shapes 
15 such as those including curves can be approximated via 
polygons that essentially flatten the curve. 

The following describes how the triangle formulae may be 
applied to three dimensional polyhedra via integrating over a 
cone. The concepts represented in FIG. 3 provide a basis for 
20 integrating over a polyhedron: 
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Lemma: Let B be a cone extended from the origin O to a 
planar (not necessarily polygonal) region F. Let h{F) be the 
distance from the origin to the plane of F. Then for any 
homogeneous quadratic polynomial / : 

5 

( 8 ) lllf(X)dv = 1 KF) jff(X)da 

B F 

Proof: As represented in FIG. 3, for any positive t (such 
as 0<r<l), the body tB is a cone extended from O to tF . For 
a small dt , the slab between tBand (t'¥dt)B is a plate defined 
10 by the region tF . The distance from O to the plane of tF is 
h{F)t , hence the slab' s thickness is h{F)dt . As an area 
integral of a homogeneous quadratic function, 
\\f{X)da^t'\\f{X)da. 

tF F 



15 



B is the integral of these slabs from 0 to 1, thus: 
\\\f{X)dv^\ \\fiX)da h{F)dt = h{F)\ \\f{X)da ydt^\h{F)\\f{X)da 



\tF 



J 



J 



A general polyhedron £ is a "signed union" of cones 
extended from O to its facets i*), where the sign is determined 
20 by the side of the plane of F^ that contains O, positive if O 
lies on the interior side. Replacing /!(/) with the 
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corresponding signed distance sh(F) ±h(f) , we have for any 



polyhedron with facets Fi,.,.,Fni 



(9) iijnX)dv^i±sh{F,)jlfiX)^^ . 




An interesting special case is when B is a simplicial 
polyhedron, i.e. when all its facets are triangles. This is 
very common in computer graphics, where hardware and software 
are highly optimized for rendering triangles. Such a 
10 polyhedron results in the following theorem and resulting 



If every facet of 5 is a triangle with 

vertices J^- J, i^-2>^-.3 ' /is a homogeneous quadratic polynomial 

15 then: 



formula: 



(10) 




Proof: Combining formulae (5) and (9) gives: 



B '=1 >=1 
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However, sh(F^)aiFi) equals three times the signed voliame of 
the cone extended from the origin to F,, and that signed volume 
equals idet(/^.„/^.2,P,^) . Thus, sh(FMP',) = \det(P,,,P,^,P.,) , hence 
formula (10) . 

A general polyhedron may be made up of the signed sum of 
the cones / pyramids. Thus, for a general polyhedron, the 
following theorem and resulting formula applies: 

Let 5 be a polyhedron with facets F,„,F„ . If / is a 
homogeneous quadratic polynomial then: 



where C, is the centroid and fl(/^)is the area of/), and 
Pi J denotes Pij-C,. 

Proof: Since sh(Fi) = -N -Ci , it follows from formulae (2) and 
(7) that: 



(11) /lJ/w^=-il;w cmfMcm 




JJJ/Wrf. = -iXW Q) /J/(^-C,) + a(is.)/(C,) , 
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and formula (7) completes the proof. 



Note that the centroid and area of can be obtained 



from: 



:i2) ti(i^.) = i|;det(i^.,., -i^„/>,, -P,M 



y«3 



Any mass properties formula can now be generated by 
defining the appropriate function / and multiplying by 
10 density. For example, }i^ithf{X)^X*X, and its polar form 
^(XpXj) = X, • ^® have from formula (4) that the moment of 
inertia of a triangle with vertices PjjPj^^j/ mass m and 
centroid C, about an axis perpendicular to the plane of T at C 
is: 

or: 



The squared distance of a point X from the axis (0,F) is 
20 (XxV)'(XxV) . This is a quadratic homogeneous function, whose 
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polar form is (X^xV)'iX2^V) . In keeping with the present 
invention, it follows from formula (5) that the moment of 
inertia of a triangle with vertices Pi^Pz^P^ and mass m about an 
axis (0,V) is: 

1=1 

Similarly, with the same squared distance function, 
multiplying formula (7) by (mass/area) density produces a 
formula for the moment of inertia of a polygonal plate about 
an axis {Q,V) , where Q is any point in the polygon's plane. 

10 Note that formula (5) produces the moment of inertia of a 

triangle about any axis through the origin, and the origin can 
be shifted to any point with a simple transformation. In 
contrast, the polygon formula will not work for an axis that 
does not intersect the polygon's plane. 

15 Again, with the same squared distance function, 

multiplying formula (10) and formula (11) by (volume/area) 
density produces formulae for the moment of inertia of a 
polyhedron about an axis (0,F) through the origin. 

The integrands for the entries of moment of inertia 

20 tensor are: 

For/j„: f{X) = y^ and the polar form isg(Jr,,X2) = >'j>'2 +^1^2 • 
For/^: f{X) = z^ and the polar form isgiX^.X^)^ z^z^-^x^x^ . 
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For/^: f(X) = x^ and the polar form isg(X^,X2) = x^X2+yly2 ' 



For/^: /(Jf) = -xy and the polar form isgiXi,X2) 



For/ 



yz 



f{X)^-yz and the polar form Is g{X^yX2) 



ForI„: f(X)--zx and the polar form is g{X^,X2) 



5 



Turning to an explanation of the present invention with 
reference to the flow diagrams of FIGS. 4 and 5^ the present 
invention can selectively combine the formulae to produce a 
desired results without integration, based on the likelihood 

10 that the most concise formula does not necessarily make the 

most efficient or most accurate algorithm. By way of example, 
the following algorithm for integrating over a polyhedron is 
likely more efficient than any straight application of 
formulae (10) or (11), described above. 

15 As represented in FIG. 4, at step 400 the integrals are 

initialized to zero. Then via step 402 {and later step 416 
which loops back to handle each facet) , for each facet with 
outward normal A^,. , if at step 404 the facet is a triangle 

(/)i,i^.2^^,3) / the process branches to step 406. 
20 Step 406 computes the area of the triangle from the 

vertices, using the above described formula a{T) = 
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j|det(P,2--'*M.-P/,3--'*/.i)i' Step 408 computes the centroid. 

With the centroid and area known, and applied to formula 
(4) described above, namely JJ/(X-C) = ^^^/(i^. , the 

5 appropriate facet integrals about centroid C,., using this 
formula are easily computed, as represented by step 410. 

Returning to step 404, if the facet is not a triangle, 
step 412 is instead executed, and computes the centroid of the 
facet using formula (13) described above, namely, the centroid 

10 of the facet, C, -j^f^^^^j-, -P^^.P^ - 

Then, at step 414, the appropriate facet integrals about Q are 
computed using formula (6) described above, namely: 

\\f{X)da = ^i;det(i^.„P,)(/-(P,_,) + /(i^) + g(P,_,.P,)) . 

With the integrals about the centroid of the facet C,. , 
15 whether computed for a triangle or otherwise, the process 
continues to FIG. 5, to compute at step 500 the facet 
integrals about the origin using formula (2a) described above, 

that is: JJ/(X)t/a = jJ/(X - QJa + a(£)/(C) . Step 502 computes the 

B B 

signed volume integrals over the facet's cone using formula 
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with N-'Cf for the signed height. 

B F 

step 506 adds the computed integrals to the accxunulative 
resulting integrals, before returning to step 416 of FIG. 4. 
As described above, step 416 of FIG. 4 returns to step 
5 402 to select the next facet with an outward normal Ni, until 
none remain, at which time the accumulative resulting 
integrals are known. It should be noted that shifting the 
polyhedron closer to the origin before executing the process 
of FIGS. 4 and 5 may further enhance accuracy. This is 

10 because when dealing with large numbers, truncation, other 

approximations and so forth are more likely to occur than when 
dealing with smaller numbers, so choosing a point as the 
origin, such as the center of mass, that will keep the numbers 
smaller, will likely be more accurate. 

15 FIG. 6 provides an example architecture 600 in which the 

present invention, such as the incorporating the algorithm of 
FIGS. 4 and 5, may be implemented. In FIG. 6, the data 602 
needed to render a graphics frame n, which may include polygon 
/ polyhedra data such as in the form of ordered facets and/or 

20 vertices 604 may be used to calculate a subsequent frame 606. 
For example, a general computation process 608 may copy the 
graphics frame 602 and modify it for the subsequent frame 606 
based on change data 610. To adjust the motion of objects to 
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compensate for collisions or other forces 602^ the general 
computation process 608 may include (or otherwise communicate 
with) the algorithm of FIGS. 4 and 5^ primarily represented by 
separate components 612-614 in FIG. 6. In general, in basic 
animation the centroid and moment of inertia is computed only 
once for a rigid body, and thus would not normally be 
recomputed for each frame. 

Note that contemporary computer graphics use large 
quantities of triangles, and thus the present invention can 
provide substantial benefits when the mass properties of an 
object or multiple objects need to be recomputed. For 
example, the moments of inertia for many smaller objects 
following the breakup of a larger object could be likewise 
computed rapidly for use in subsequent frames using the 
present invention. As another example, an image of an ice 
skater spinning faster by moving in his or her arms (which 
changes the moment of inertia) could be rapidly calculated for 
each frame to provide more realistic rendering. 

As can be seen from the foregoing detailed description, 
there is provided a method and system to compute the mass 
properties of polygons and polyhedra, including by which the 
moment of inertia of polygons and polyhedra may be efficiently 
and accurately determined. The method and system also extend 
to any quadratic homogeneous polynomial. The present 
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invention may be used as a convenient and straightforward 
mechanism for determining mass properties for graphics 
processing, as well as in other technical applications, and 
thus provide significant advantages and benefits needed in 
5 contemporary computing. 

While the invention is susceptible to various 
modifications and alternative constructions, certain 
illustrated embodiments thereof are shown in the drawings and 
have been described above in detail. It should be understood, 
10 however, that there is no intention to limit the invention to 
the specific forms disclosed, but on the contrary, the 
intention is to cover all modifications, alternative 
constructions, and equivalents falling within the spirit and 
scope of the invention. 
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